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Abstract. We present the results of a search for Be stars 
in six fields centered on the young clusters NGC 330 and 
NGC 346 in the SMC, and NGC 1818, NGC 1948, NGC 
2004 and NGC 2100 in the LMC. Be stars were identi- 
fied by differencing R band and narrow-band Ha CCD 
images. Our comparatively large images provide substan- 
tial Be star populations both within the clusters and in 
their surrounding fields. Magnitudes, positions and finding 
charts are given for the 224 Be stars found. The fraction 
of Be stars to normal B stars within each cluster is found 
to vary significantly although the average ratio is similar 
to the average Be to B star ratio found in the Galaxy. In 
some clusters, the Be star population is weighted to mag- 
nitudes near the main sequence turn-off. The Be stars are 
redder in V—I than normal main-sequence stars of similar 
magnitude and the redness increases with increasing Ha 
emission strength. 

Key words: Stars: emission-line - Stars: Be - Galaxies: 
clusters - Galaxies: Magellanic Clouds 



1. Introduction 

In this paper we present the results of a photometric sur- 
vey for Be stars within and around six young clusters in 
the Magellanic Clouds. Such studies have been carried out 
on Galactic clusters but they are limited there by the 
relative paucity of Be stars. The young populous Mag- 
ellanic Cloud clusters examined in this paper present us 
with large samples of Be stars in a range of metallicity 
environments and the opportunity to further examine the 
mechanism for formation of Be stars. 

Searches for Be st ars in the Magellanic Clouds have 
been made by Feast ( 1972 ) using objective prism tech- 
niques, and the m ore r ecently by Grebel et al. ( 1992| ), 
Bessell and Wood ( 1992 ) (who gave a brief description of 
s ome of the results presented in this paper) and Grebel 
(1997) using CCD-based imaging photometry. These lat- 
ter searches involved imaging through an Ha filter so that 
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Be stars, which show strong Ha emission, stand out in 
comparison to normal stars. Photometric techniques such 
as these are very efficient methods of identifying Be stars, 
in particular within dense clusters where spectroscopy is 
difficult. 

The main purpose of the present study is to provide 
a sample of Be stars for follow-up spectroscopy and to 
see if there are differences in Be star populations between 
different clusters and between clusters and the local field 
population. The results of the follow-up spectroscopy will 
be presented in a subsequent paper. Because the present 
study is ground-based, its spatial resolution is seeing lim- 
ited, which means that we are unable to search for Be stars 
in the cluster cores. In another paper, we will present the 
results of a search for Be stars in the cluster cores using 
HST imaging. The present paper provides the comparison 
sample of Be stars in the local field and the outer cluster 
zones. The complete data set of ground-based and HST 
images will allow an accurate comparison of field and clus- 
ter Be stars properties. In this paper, however, we make 
some initial comparisons with the more limited data set 
described herein. 

2. Observations and Data Reduction 

The photometric search method we used to detect Be stars 
involved taking a CCD image through a narrow-band (15 
A) Ha filter and comparing this image with a similar im- 
age obtained in the Cousins R band. Stars with strong Ha 
emission should appear brighter in the Ha image than in 
the R image. Given that a typical Be star has Ha emission 
with a full width at half maximum of ^7 A and peak Ha 
flux 5 times the local continuum flux, this sort of search 
with a 15 A filter is readily capable of detecting Be stars. 

The observations were obtained using direct CCD 
imaging on the 1-m telescope at Siding Spring Observa- 
tory from September to November 1991. All Ha images 
were 900 second exposures. R images were obtained im- 
mediately before or after each Ha image in order to avoid 
the possibility of stellar variability affecting our measure 
of emission strength. A Tektronix 1024x1024 CCD with 
pixels of scale 0.6" was used for these observations, giving 
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Table 1. The log of observations presented in the current 
paper. 



Cluster 


Date(1991) 


Filter 


Exposure Time 


NGC 330 


Sept. 04 


V 


300s 






I 


150s 




Oct. 31 


R 


300s 






Ha 


900s 


NGC 346 


Nov. 29 


V 


150s 






I 


100s 






R 


150s 






Ha 


900s 


NGC 1818 


Sept. 04 


V 


80s 






I 


150s 




Oct. 31 


R 


300s 






Ha 


900s 


NGC 1948 


Oct. 30 


V 


120s 






I 


120s 






R 


200s 






Ha 


900s 


NGC 2004 


Nov. 29 


V 


20s 






V 


150s 






I 


15s 






1 


100s 






R 


15s 






R 


100s 






Ha 


900s 


NGC 2100 


Nov. 29 


V 


120s 






I 


100s 






R 


100s 






Ha 


900s 



us an imaged area of typically 10'xlO' (compared with 
5.7'x5.7' in Grebel 1997). V and / images were obtained 



at the same time for some fields, and in some cases on 
other nights using a GEC 770x1152 CCD. A log of the 
observations is presented in Table [l| 

The CCD images were processed with IRAF and the 
photometry of the fields was done using the DoPHOT pho- 



tometry package (Mateo and Schechter 1989). Standard 
magnitudes in V and / were derived from the standards 



used in Sebo and Wood (1994) and from magnitudes given 
by Walker (1996). The R and Ha magnitudes were not 
standardised and R— Ha colours have an arbitrary zero- 
point. 

3. Selection of Be and other Ha emitting stars 

In order to detect stars with Ha emission, plots were made 
of the colour index R— Ha against V—I for all the stars 
which were common to the V, I, R and Ha frames. Figs, [j] 
(NGC 330), | (NGC 346), | (NGC 1818), | (NGC 1948), 
| (NGC 2004) and | (NGC 2100) show the (i?-Ha, V-I) 
diagrams produced. As noted in Section [| the R— Ha 
zero-point is arbitrary in these plots, although the zero- 
point was chosen so that normal main-sequence stars had 
i?-Ha w 0. 
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Fig. 1. The (V, V-I) and (R-Kct, V-I) diagrams for 
the cluster NGC 330 (panels a and c) and its surrounding 
field (panels b and d). The usable cluster area lies between 
radii 15" < r < 1'.6 while the field is defined as r > 1'.6. 
Be stars are identified as the stars above the dashed line 
in the (R— Ha, V—I) diagrams and they are shown on the 
(V, V—I) diagrams as filled circles. 
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Fig. 2. The same as Fig. 0, but for NGC 346. The usable 
cluster area lies between radii 15" < r < 2'. 2 while the 
field is defined as r > 2'. 2. 



The stars with strong Ha emission were selected from 
the above figures. Using Fig. [l]c as a typical case, we note 
that most main-sequence stars lie in a tight clump around 
R— Ha=0.0 and V— 1=— Q.15, while nearly all cooler stars 
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Fig. 3. The same as Fig. [j], but for NGC 1818. The usable Fig. 5. The same as Fig. [I], but for NGC 2004. The usable 
cluster area lies between radii 15" < r < 2'. 2 while the cluster area lies between radii 15" < r < 2'.0 while the 



field is defined as r > 2'. 2. 



field is defined as r > 2'.0. 
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Fig. 4. The same as Fig. @, but for NGC 1948. The usable Fig. 6. The same as Fig. Q, but for NGC 2100. The usable 
cluster area lies between radii 15" < r < 2'. 7 while the cluster area lies between radii 15" < r < 1'.4 while the 



field is defined as r > 2'. 7. 



field is defined as r > l'A. 



form an almost horizontal band. These two groups of stars 
do not have detectable Ha emission. However, there is 
another group of stars extending up and slightly to the 
red of the main-sequence clump: these stars clearly show 
significant Ha emission and they are Be stars. Occasional 
cooler (giant) stars also exhibit strong Ha emission. 

The dispersion in R— Ha of normal MS stars about the 
centre of the clump is determined largely by the photo- 
metric errors in the narrow-band Ha magnitude, and this 



uncertainty will depend on the magnitude of the star. We 
limit our Be star selection process to stars with V<17.5, 
which corresponds approximately to an error in R— Ha 
colour of ^0.2 magnitudes. In order to select a sample of 
Be stars, we draw a line in each of the R— Ha diagrams 
parallel to the sequence of non-emission line stars and at a 
distance of ^0.4 magnitudes above it, which gives a high 
probability that we exclude all members of the clump of 
normal main-sequence stars. 
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Our observations were made under seeing conditions 
of 1.5-2", and as a consequence the central regions of the 
clusters suffer from crowding. We therefore exclude the 
innermost 15" (in radius) of each cluster from the Be star 
search. On each image taken, we select a radius at which 
the stellar density drops to a value indistinguishable from 
that of the surrounding field. Stars within this radius are 
assumed to belong to the cluster, and those outside to the 
field. The radius selected for each cluster is given in the 
captions of Figs. |l| to [| 

The positions, magnitudes and colours of the Be stars 
selected as described above are listed in Tables 2 (NGC 
330), 3 (NGC 346), 4 (NGC1818), 5 (NGC 1948), 6 (NGC 
2004) and 7 (NGC 2100). Tables 2-7 are available in elec- 
tronic form from the CDS. The astronomical coordinates 
in these tables were derived from the secondary astromet- 



.2. NGC 1818 
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ric standards of Tucholke et al. (1996). Finding charts for 
the Be stars are given in Figs. 
|H and|l7| (NGC 346), [jj and ~ 



and h5| (NGC 330), 



(NGC 1818), EOl (NGC 



1948) , fc2_y and E2| (NGC 2004), and and 24 (NGC 2100) 



4. Comparison with past surveys for Be stars 

We now compare the results of our study to those of pre- 
vious Be star surveys within NGC 330 and NGC 1818 
in order to get some idea of our detection efficiency. Al- 
though there have been a few reports of detecti ons of Be 
stars in other clusters (e.g. Kjeldsen and Baade 1994 ), it 
is only for NGC 330 and NGC 1818 that any information 
is available regarding the identity of the Be stars. 

4.1. NGC 330 

Grebel et al. ( |1992|) have searched the region close to 
NGC 330 for Be stars using narrow-band Ha photome- 
try. Within 50" of the cluster centre they find 29 Be stars 
while we find 32. However, the Be star samples found in 
the two respective searches are not completely identical. 
Of the 18 Be stars for which Grebel et al. provid e identi- 
fications (they list names from Robertson 1974 and Arp 



195E), we find 16. We suspect that variability of the Be 



phenomenon may account for the differences in the two Be 
star samples (see Section 12). Overall, however, the two 
searches find similar numbers of stars. Our search area is 
considerably larger than that of Grebel et al. which will 
allow us to compare the properties of reasonable samples 
of both field and cluster Be stars. 

NGC 330 and its surroundings have been searched ex- 
tensively for variable stars by Sebo and Wood (1994) and 
Balona (1992). A number of the Be stars and one of the 
Ha emitting red giants found here are known variables. 
The variable star identifications are listed in Table 2. A 
spectrum of the v ariable red giant was obtained by Sebo 
and Wood ( 1994 ) and it does indeed show a broad Ha 
emission line. 



Grebel ( 1997 ) presents the results of a Be star search 
within a 5.7' x 5.7' field around NGC 1818. Outside a ra- 
dius of 15" from the cluster centre, a total of 48 Be stars 



were found. The study of Grebel ( 1997 ) was made under 
excellent seeing conditions and has a limiting magnitude 
of V ~ 20. Grebel's selection criterion required Be stars 
to be ^0.2 magnitudes above the main-sequence clump in 
the (R— Ha, V— I) diagram, compared to our ~0.35 mag- 
nitudes (see Grebel's Fig. 3 with our Fig. ||b). However, 
given our Ha filter bandpass of 15 A and Grebel's band- 
pass of 33 A, the sensitivities of the two surveys should be 
similar, at least for the brighter stars where photon noise 
in the Ha filter is ^0.35 magnitudes in our case, or 0.2 
magnitudes in Grebel's case. 

Grebel fll997| ) finds 21 Be stars with V < 17.5 (the 
magnitude limit of our study) , of which we have identified 
13 in common. Of the 8 unidentified stars, 7 are near our 
magnitude limit and have 17.1 < V < 17.5: furthermore, 
most lie in crowded regions. More interestingly, we de- 
tected 6 Be stars that have no counterparts in the survey 
of Grebel. This is perhaps surprising given the fainter mag- 
nitude limit and better seeing of the Grebel survey. These 
stars, which are all comparatively bright and lie in un- 
crowded regions, are unlikely to be misidentifications. We 
thus seem to have strong evidence for the episodic nature 
of the Be phenomenon in Magellanic Cloud stars: such be- 
haviour is frequently seen in galactic Be stars (Hanuschik 
et al. 1993). It is also possible that some of the Be stars 
found by Grebel, but which we failed to detect, were not 
exhibiting Ha emission at the time of our observation. 

In summary, both the present survey and the surveys 
of Grebel et al. (1992) and Grebel (1997) find similar num- 
bers of Be stars, although our survey may be incomplete 
near the magnitude limit of V = 17.5 in crowded regions. 
Because of the episodic nature of the Be phenomenon, 
single-epoch surveys may miss a significant fraction of 
the total Be star population (~20% appears to have been 
missed in the case of NGC 1818). 

5. Photometric properties of the Be star 
population 

There are a number of interesting correlations between the 
photometric indices of Be stars observed here which may 
shed some light on the nature of the Be star phenomenon. 
Starting with Figs. |j]to^| we see that within the clusters, 
as well as in the field, Be stars exist over a wide range of 
luminosities. Since the Be stars in the clusters are all of 
essentially the same age (field contamination of the cluster 
Be star populations should be <10% based on the small 
cluster area), the wide range in Be star luminosities means 
that the Be phenomenon can occur at any time throughout 
the main-sequence lifetime, and need not be confined to 
particular evolutionary phases such as the core contraction 
phase at the end of the main-sequence. 
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Another interesting correlation is seen in the (R— Ha, 
V—I) plots (Figs. 0c to ^|c). Here we see a significant trend 
for the strongest emitters, as measured by the R— Ha 
colour, to be those with the reddest V—I. This effect was 



also noted by Grebel (1997). The fact that Be stars are 
redder in V—I than normal B stars is also obvious from 
Figs. |l|a,b to ||a,b where, at a given V, the Be stars are 
clearly redder than the normal main-sequence stars. An 
increase in redness of Be stars with increased strength 
of Ha emiss ion h as also been noted in the infrared (e.g. 
Dachs et al. |l98S| ). 

It is widely accepted that the Be phenomenon is asso- 
ciated with rapid rotation of the stellar photosphere and 
the presence of a circumstellar disk of comparatively cool 
material which gives rise to the observed Ha emission. 
The observed increase in V—I of Be stars with increased 
strength of Ha emission is almost certainly the result of 
one or more of four factors: a real reduction in stellar 
effective temperature due perhaps to rotation, a changed 
spectral energy distribution due to rotational distortion of 
the stellar atmosphere, continuum and line emission from 
the disk, or absorptive reddening caused by circumstellar 
matter. The first two factors result from changes in the 
structure of the star, while the latter two involve circum- 
stellar matter only. The viewing angle will also effect the 
colours in the non-spherical geometry. 

In this regard our photometric colours do not consti- 
tute a conclusive test of the mechanism responsible for the 
apparent reddening of Be stars. Fig. shows plots of B—V 
and V — R against V—I for Be stars and normal B stars 
within and around the cluster. The B—V values come from 
Balona and Jerzykiewicz ( |1993[ ). The zero point of V—R 
is arbitrary (since the R magnitude has not been zero- 
pointed) but has been adjusted to be close to the value 
expected for the normal hot B stars. 

In the (V—R, V—I) plot, the Be stars fall on a se- 
quence that is essentially indistinguishable from the se- 
quence followed by non-Ha emitting stars. This overlap is 
coincidental and is the consequence of the spectral distri- 
bution of the continuum emission from the circumstellar 
within the R and / bands. 



envelope (Kaiser 1989 ) within the R and / bands. The 
position of the Be stars in the (B—V, V—I) in Fig. is 
suggestive of a systematic shift to the right of the sequence 
of normal main-sequence stars, shown by the continuous 
line. In this case significant emission due to Paschen con- 
tinuum emission from the circumstellar envelope is located 
in the / band alone. Hence the Be stars appear shifted to 
the right relative to the non-Be stars. 

Also shown in Fig. are the reddening lines from 
Taylor fll 98^ ). On the (B—V, V-I) plot, it is clear that 
the Be stars do not follow the reddening line, suggesting 
that the red colours of the Be stars are not due to circum- 
stellar dust absorption. 

The effect on B—V colour of rotational distortion of 
the at mosphe re of a Be star has been examined by Collins 
et al. ( 1991 ). However, even in the case of a star ro- 
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Fig. 7. The (V-R, V-I) and (B-V, V-I) diagrams 
for stars within and around the cluster NGC 2004, and 
brighter than V=16.5. Non-emission B stars are shown as 
dots and Be stars are shown as solid squares. The B—V 



colours are taken from Balona and Jerzykiewicz (1993). 
Due to the relatively small field studied by Balona and 
Jerzykiewicz, there are fewer stars in the (B—V, V—I) 
diagram. The solid line in the (B—V, V—I) diagram is 
the theoretical main-sequence (logg — 4.5) from Bessell 
et al. ( 199Sp . The arrows are interstellar reddening vec- 
tors from Taylor ( 1986| ), corresponding to E(B— V) = 0.1. 



tating near breakup, the maximum predicted change in 
B—V colour (from the colour of a no n-rota ting B star) is 
only ^0.05 magnitudes (Zorec & Briot 1997 ), compared to 
changes of up to 0.15 magnitudes seen in Fig. 0. It there- 
fore seems that this effect is not a major contributor to 
the reddening of Be stars. 

The final possibility is that the red colours of the Be 
stars compared to the colours of non-Be stars at the same 
luminosity (as seen in Figs. 0a to ^a) are due to a real 
change in the photospheric flux distribution through the 
BVRI bands caused by a change in stellar effective tem- 
perature. However, since the bulk of the flux of these stars 
is emitted shortward of B, it remains to be seen whether 
the Be stars really do have lower effective temperatures 
than B stars of similar luminosity. This question will be 
examined in forthcoming papers which will give spectral 
types of the stars, and HST ultraviolet magnitudes from 
which T c g values can be accurately compared between the 
Be and normal B stars. 
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6. Comparison of cluster and field Be star 
populations 

The Magellanic Clouds offer us an opportunity to examine 
the held population of Be stars in different localities in 
an unbiased way - a task which is virtually impossible 
from within our own galaxy due to depth and reddening 
effects. For the purposes of our survey we can regard each 
CCD image as essentially a volume and magnitude limited 
sample of the population. 

We have defined the cluster extent by a direct exami- 
nation of the radial variation in surface density of bright 
blue stars (V < 17 and V— I < 0.6), as described in Sec- 
tion H Some stars co-eval with the cluster stars proba- 
bly exist outside o ur ad opted radii. In particular, we note 
that Elson et al. ( 1987 ) have examined the surface den- 
sity of stars around a sample of young LMC clusters and 
suggested that some of these clusters are tidally unbound 
with up to 50% of their total mass in an unbound halo (al- 
though the halo is of very low surface density compared 
to the background). 

The CMDs for both the cluster and field of NGC 330, 
NGC 346, NGC 1948, NGC 2004 and NGC 2100 are very 
similar (see Figs. [TJd to J^d). This indicates that the field 
populations contain stars of age similar to that of the clus- 
ter. We have investigated the fraction of stars along the 
main-sequence that arc Be stars by binning both emission- 
line and non-emission-line stars with V—I < 0.5 into 1/2 
magnitude bins in V. Figs. || to |l3| show our results. In Ta- 
ble ^, we give the numbers of Be stars and main-sequence 
stars with V < 17 for each cluster and its surrounding 
field: we also give the ratio of Be to main-sequence stars 
(with a l-a error). 

Here we have made the assumption that comparison 
of the number of Be stars to B stars within a certain lu- 
minosity range is a valid comparison. Th is ass umption is 
challenged by the study of Zorec & Briot ( 1997 ) who claim 
that the location of Be stars within the CMD is the result 
of a visual luminosity excess. The V excess described by 
Zorec & Briot rises from 0.0 at B9 to -0.5 at B0 in a linear 
manner. Here we do not adjust for V excess as we consider 
the size of the associated uncertainties in the correction 
similar to the size of the correction itself. 

The fraction of Be stars in the field populations lies 
in the range 0.10 to 0.27. The mean value is a similar 
to the fraction of Be stars in the Galactic field which is 



0.17 (Zorec and Briot 1997). The fractions of Be stars in 
the clusters is quite variable. The clusters NGC 330 (see 
also Grebel et al. |1992[ ) and NGC 2100 are particularly 
rich in Be stars, NGC 1818 has a moderate richness, while 
NGC 346, NGC 1948 and NGC 2004 have low Be star 



fractions. In the Galaxy, Zorec and Briot (1997) noted 
that the Be star fraction reached a maximum of 0.34 at 
spectral type Bl. The fraction of Be stars in NGC 330 
(0.34±0.08) is perfectly consistent with this value. In a 
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Fig. 8. a. Histograms of the number of Be stars (shaded) 
and main-sequence stars (V—I < 0.5, including Be stars) 
in half-magnitude bins down the main-sequence of NGC 
330. b. Same as a but for the field around NGC 330. c. The 
ratio of Be stars to main-sequence stars in half-magnitude 
bins down the main-sequence of NGC 330. d. Same as c 
but for the field around NGC 330. 



forthcoming paper, the correlation of Be star ratio with 
spectral type will be examined. 

In order to compare the Be star fraction near the main- 
sequence turnoff with the Be fraction for less evolved stars, 
we have divided the stars into two groups: the brighter 
stars with V < 16 and the fainter stars with 17 < V < 16, 
all stars having V—I < 0.6. The Be star fractions for these 
groups are given for the clusters and their surrounding 
fields in Table |^. From these ratios, and looking at Figs. |^ 
to |l^, it seems that for most of the clusters and fields there 
is no statistically significant difference between the Be star 
fractions near the main-sequence turnoff and further down 
the main-sequence. There are a few distinct exceptions to 
this rule: the clusters NGC 330 and NGC 2004 show a 
significantly higher fractions of Be stars near the main- 
sequence turnoff than at fainter magnitudes, as does the 
field of NGC 1818 (where the four brightest main-sequence 
stars are all Be stars). 

7. Discussion 

We now examine what our observations imply about the 
Be mechanism. In the literature, binarity, pulsation and 
rotation have been suggested as the mechanism for the 
establishment and maintainence of the circumstcllar ma- 
terial seen in Be stars. 
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cluster stars 


field 


stars 




V < 17 


16 < V < 17 


V < 16 


V < 17 


16 < V < 17 


V < 16 


Cluster 


Nms 


N Bc 


Nbo/Nms 


Nbc/Nms 


Nbc/Nms 


Nms 


N Bc 


Nbo/Nms 


Nbc/Nms 


Nbc/Nms 


NGC 330 


79 


27 


0.34±0.08 


0.25±0.08 


0.48±0.15 


128 


27 


0.21±0.04 


0.12±0.04 


0.19±0.08 


NGC 346 


78 


11 


0.14±0.05 


0.19±0.08 


0.11±0.05 


127 


23 


0.18±0.04 


0.19±0.06 


0.17±0.05 


NGC 1818 


92 


19 


0.21±0.05 


0.16±0.06 


0.27±0.10 


38 


10 
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Fig. 9. Same as Fig. § but for NGC 346. 



Fig. 10. Same as Fig. | but for NGC 1818. 



Harmanec and Kfiz ( 1976 ) have suggested that Be 
stars comprise a population of interacting binaries. In this 
model, elaborated upon by Pols et al. ( 1991 ), emission 
arises from material in accretion onto the secondary com- 
ponent. The difficulty that remains with this model is that 
the expected number of close interacting binaries is too 
few to account for the observed population of Be stars 



(VanBever and Vanbeveren 1997). Our observations, par- 
ticularly the high Be star ratio in NGC 330, do not favour 
the interacting binary hypothesis as the population of 
stars required to be in mass exchanging binary systems 
near the MS turnoff (Abt |1987| ) is too high. 

The non-radial pulsator model for Be activity (Vogt 
and Pcnrod 1983; Baade 1987) has gained support from 
the recent work of Dziembowski et al. (1993). In their 



models, Dziembowski et al. have shown that an instability 
domain exists on the main-sequence that almost fills the 
gap between S Set and (3 Cep variables. These slowly pul- 
sating B-type stars span a range of mass from 3 to 6 M . 
The B e star (i Cen has been intensively studied by Rivinius 
et al. ( 1998 ). From their analysis of photospheric radial ve- 



locity variations they have discerned six discrete periods 
between 0.27-0.51 days. Rivinius et al. suggest that when 
the majority of these modes are in constructive interfer- 
ence, expulsion of material and a line emission outburst 
occurs. 

There are problems with associating the Magellanic 
Cloud Be stars observed here with the pulsation instabil- 
ity mechanism described above. Firstly, in conditions of 
reduced metallicity such as those within the Magellanic 
Clouds, the extent of the pulsation instability domain is 
reduced (Dziembowski et al. 1993[) and f3 Cep instability is 
seemingly extinguished (Balona 1992). Most importantly, 
the Be stars observed here are by no means confined to 
a range 556 M Q . Indeed our observations show that the 
largest proportion of Be stars occurs at ~ 14-15 M Q . 

The most widely accepted explanation for the Be phe- 
nomenon is that it is a result of the rapid rotation of the 
B star. This explanation was originally put forward by 
Struve (1931). In it s mo dern incarnation, due to Bjork- 



man and Cassinelli (11993D , the stellar wind from a rapidly 
rotating star is concentrated and confined to a disk around 
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Fig. 13. Same as Fig. | but for NGC 2100. 
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Fig. 12. Same as Fig. | but for NGC 2004. 



the equatorial plane of the star. The model of Bjorkman 
and Cassinelli predicts a maximum likelihood for disk for- 
mation near spectral type B2 (where the stellar angular 
velocity required for disk formation is 55% that of the crit- 
ical breakup velocity). Observations of Galactic Be stars 
find that the highest Be star fraction does indeed occur 
at ea rly spectral types (Bl according to Zorec and Briot 
1997 ). Within our sample spectral types are available for 



only five of our flagged Be candidates. We will present 
further spectral types in a forthcoming paper. 

Assuming the Be phenomenon is related to the forma- 
tion of an equatorial disk caused by rapid rotation, we 
would expect that clusters which form from gas that has 
a high angular momentum content would contain a high 
fraction of Be stars. Indeed, the evidence is suggestive that 
Galactic clusters whose stars present a higher than aver- 
age rotational velocity have larger populations of Be stars 
(Mermilliod 1982). In the Magellanic Clouds, the cluster 
NGC 330 has a noticeably ellipsoidal shape which may 
suggest that the cluster as a whole formed from material 
of higher angular momentum and that, as a consequence, 
the cluster stars have a rapid rotation rate. This would 
certainly be consistent with the high Be star content of 
NGC 330. In a future paper we will examine the rotation 
velocities of samples of stars within and around the clus- 
ters NGC 330 and NGC 2004 to see if the overall Be star 
fraction is related to average cluster star rotation velocity. 

The tendency for there to be a higher Be star frac- 
tion near the main-sequence turnoff in some clusters also 
appears consistent with the rapid rotation hypothesis for 
Be star formation. The evolution of rapidly rotating stars 
has been examined by Endal and Sofia (1979) and Endal 
(1982). These models explicitly follow the radial exchange 
of angular momentum with evolution. 

In the models of Endal and Sofia, stars commencing 
their main-sequence lives with relatively slow angular ve- 
locity remain slow rotators throughout the whole of the 
main-sequence evolutionary phase (although there is a 
marked increase in angular velocity of all stars during the 
core contraction phase associated with exhaustion of hy- 
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drogen in the core, this evolutionary phase is far too short 
to account for the observed proportion of Be stars). Thus 
slow rotators will never evolve into Be stars. 

In contrast to the slow rotators, the models of En- 
dal and Sofia show that stars commencing their lives with 
an angular velocity greater than 56-76% of the critical 
breakup velocity lo c spin up to the critical velocity over a 
moderate fraction of the main-sequence lifetime. In a clus- 
ter where some fraction of the stars is formed with angular 
velocity greater than ~50% of u> c> those stars nearer the 
main-sequence turnoff are more likely to be Be stars since 
they will have been more spun up by evolution than the 
less luminous stars that have not evolved as far through 
the main-sequence phase. This process could explain the 
concentration of Be stars to the main-sequence tip seen 
in NGC 330, NGC 2004 and NGC 1818. We add the cau- 



tionary note that Endal and Sofia Ql979| ) and Endal ( |1982| ) 
made a detailed study of only a single star which had a 
mass of 5 M©, much less than the turnoff masses of ^15 
Mq which is appropriate for the clusters observed here. 
For comparison with Magellanic Cloud data, models such 
as those of Endal and Sofia need to be made at higher 
masses and lower metallicities. 

Other factors which might influence the Be star frac- 
tion are metallicity and age. However, neither of these 
factors seem to influence the Be star fractions in our 
data. To examine abundance dependence, we note that 
the two SMC clusters are probably more metal poor than 
the LMC clusters by a factor of about two (Russell and 
Bessell pJ89| ). However, in the SMC, NGC 330 has a high 
Be star content while NGC 346 does not. At the higher 
abundance of the LMC, clusters with a range of Be star 
content are present, and the range is similar to the range 
in the SMC in spite of the abundance difference between 
SMC and LMC. Turning to age, we note that the clus- 
ters studied here have ages of ~l-3xl0 7 years, although 
these ages are uncertain by factors of 2-3 because of uncer- 
tainties in the amount of convective core overshoot which 
occurs in real stars. For comparison of Be star fractions 
with age, we take from the literature the following ages: 
NGC 330, 20 Myr (Chiosi et al. |1995| ); NGC 346, 1 3 My r 
(Massey et al. |1989|) ; NGC 1818, 30 My r (Wi ll et al. |1995| ); 
NGC 1948, 15±5 Myr (Vallen ari et al. |1993| ); NGC 2004, 
20 Myr (Caloi a nd Ca ssatella |1995| ); NGC 2100, 15 Myr 
(Cassatella et al. 1996). There does not appear to be any 
correlation of Be star fraction with age: the cluster NGC 
330 with the highest Be star fraction has an intermediate 
age and the oldest cluster (NGC 1818) has an intermedi- 
ate Be star fraction. This result contrasts with the results 



of Mermilliod (1982) who noted a trend of decreasing Be 
fraction with increasing age amongst young Galactic clus- 
ters. 



8. Summary 

We have presented the results of a photometric survey for 
Be stars in six fields centred on NGC 330 and NGC 346 
in the SMC, and NGC 1818, NGC 1948, NGC 2004 and 
NGC 2100 in the LMC. In all six of these young populous 
clusters we find large numbers of Be stars. 

The Be stars are shown to be redder in V— I than 
main-sequence stars of similar My. Furthermore, the dis- 
placement to the red increases with increasing Ha emis- 
sion line strength. Some suggested causes for the redness of 
Be stars were examined using the widely accepted model 
that Be stars are main-sequence stars with a surrounding 
disk. Our photometric observations suggest that the red- 
ness of the stars is due to emission from the circumstellar 
envelope, in line with the results of detailed spectropho- 
tometric studies based upon galactic samples. 

The fraction of main-sequence stars that are Be stars 
varies significantly between the various clusters from 
~0. 10-0.34. The average value is close to the average frac- 
tion of ~0.17 seen in the Galaxy, and the maximum value 
is also similar to the maximum Galactic value, which oc- 
curs at spectral type Bl. It was noted that the cluster 
with the maximum Be star fraction, NGC 330, is dis- 
tinctly elliptical, suggesting that the stars in it formed 
out of rapidly rotating material and that there is a con- 
nection between rapid rotation and the Be phenomenon. 
No connection between Be star fraction and age or metal- 
licity was evident. In some clusters, the Be stars tend to 
be concentrated toward the main-sequence turnoff. This 
is consistent with evolutionary models which predict that 
stars which form with rotation velocities more than ^50% 
of the breakup value spin up as they evolve, reaching the 
breakup velocity over a moderate fraction of the main- 
sequence lifetime. 
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Fig. 14. The finding chart (R image) for the Be stars around NGC 330. The full field is 10'xlO' while the inner 
rectangle is 7.'3 x 7.'9 and is the part of the field on which the search for stars with Ha was performed. North is to the 
right and east is at the top. 
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Fig. 15. Same as Fig. |lj but for stars near the core of NGC 330. 
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Fig. 16. The finding chart for the Be stars around NGC 346. The field size is 10' square and the orientation is north 
to the right, east at the top. 
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Fig. 17. Same as Fig. [ll] but for stars near the core of NGC 346. 
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Fig. 18. The finding chart for the Be stars around NGC 1818. The full field is 10'xlO' whilst the inner rectangle in 
which the search for stars with Ha emission was performed is 7.'6 x 9.'2. North is to the right and east is at the top. 
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Fig. 19. Same as Fig. [T^ but for stars near the core of NGC 1818. 
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Fig. 20. The finding chart for the Be stars around NGC 1948. The full field is 10'xlO' whilst the inner rectangle in 
which the search for stars with Ha emission was performed is 7.'6 x 9.'2. North is to the right and east is at the top. 
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Fig. 21. The finding chart for the Be stars around NGC 2004. The field size is 10'xlO'. North is to the right and east 
is at the top. 
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Fig. 22. Same as Fig. |l] but for stars near the core of NGC 2004. 
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Fig. 23. The finding chart for the Be stars around NGC 2100. The field size is 10'xlO'. North is to the right and east 
is at the top. 
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Fig. 24. Same as Fig. |23| but for stars near the core of NGC 2100. 
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